function res = run_iodine_reaction()

% initial concentration
I = .5     ;     % V(1)
S2O8 = .5;        % V(2)
I2 = 0;         % V(3)
SO4 = 0;        % V(4)
S2O3 = .25;       % V(5)
S4O6 = 0        % V(6)
endTime=10; 

[t, V] = ode45(@iodine_reaction, [0,endTime], [I, S2O8, I2, SO4, S2O3, S4O6]);

plotter(t,[V(:,5) V(:,3)],...
    'Concentration over time for S203 and Blue Complex', 'Time (s)', 'Concentration',...
    ['S203        '; 'Blue Complex'],[ [1 0 0];[ 0 0 1]])

end
